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Abstract 



Reliable values of quark and lepton masses are important for model building 
at a fundamental energy scale, such as the Fermi scale M z ~ 91.2 GeV and 
the would-be GUT scale A GUT ~ 2x 10 16 GeV. Using the latest data given by 
the Particle Data Group, we update the running quark and charged-lepton 
masses at a number of interesting energy scales below and above M z . In 
particular, we take into account the possible new physics scale (/i ~ 1 TeV) 
to be explored by the LHC and the typical seesaw scales (/i ~ 10 9 GeV and 
/i ~ 10 12 GeV) which might be relevant to the generation of neutrino masses. 
For illustration, the running masses of three light Majorana neutrinos are also 
calculated. Our up-to-date tables of running fermion masses are expected to 
be very useful for the study of flavor dynamics at various energy scales. 
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I. INTRODUCTION 



In the standard model (SM) of electroweak interactions [1], it is the Higgs mechanism 
that provides a self-consistent framework to generate masses for gauge bosons and charged 
fermions. Three neutrinos are exactly massless as a straightforward consequence of the 
symmetry structure of the SM. But this framework itself can neither predict the values of 
six quark masses and three charged-lepton masses nor interpret the observed hierarchy of 
their spectra [2]. On the other hand, current neutrino oscillation experiments [3-6] indicate 
that neutrinos are actually massive and lepton flavors are mixed. Hence one has to extend 
the SM in order to gain an insight into the dynamics of fermion mass generation and flavor 
mixing. Possible new physics beyond the SM is expected to appear far above the Fermi scale 
M z 91.2 GeV which essentially characterizes the scale of electroweak symmetry breaking. 
For instance, the typical scale of grand unified theories (GUTs) for strong and electroweak 
interactions [7] would be A GUT ~ 2 x 10 16 GeV, while /i ~ 1 TeV is the possible scale of new 
physics necessary to stabilize the mass of the Higgs boson and offer a natural explanation 
of electroweak symmetry breaking [8]. Whenever a specific mass model of leptons and 
(or) quarks is built at a certain high energy scale, one has to make use of the technique of 
renormalization group equations (RGEs) [9] to bridge the gap between the model predictions 
(at fi ^> M z ) and the experimental data (at \i < M z ). Therefore, reliable values of running 
quark and lepton masses are very important for building new physics models and testing 
their viability. 

A systematic estimate of running quark masses at various energy scales was made by 
Fusaoka and Koide [10] in 1998. Their numerical results have proved to be very useful for 
people working on both hadronic physics at relatively low energies and electroweak physics 
and GUTs at superhigh energies. However, a nontrivial update of this work is greatly 
desirable today because a lot of changes in particle physics have happened since 1998. On 
the one hand, some new and more accurate data on fermion masses at low scales have been 
accumulated [11]. On the other hand, intensive interest has been paid to some new physics 
scales such as /i ~ 1 TeV to be explored by the Large Hadron Collider (LHC) and /x ~ 10 12 
GeV which might be relevant to the seesaw mechanism of neutrino mass generation [12]. 
These energy scales were not considered in Ref. [10], nor were the running masses of three 
light neutrinos. Thus we are well motivated to do a new analysis of running fermion masses 
at a variety of fundamental or interesting energy scales. 

Let us remark that the present paper is different from Ref. [10] and other previous works 
in the following aspects: 

• The running masses of three light Majorana neutrinos are calculated. A global analysis 
of current neutrino oscillation data yields strong constraints on two neutrino mass- 
squared differences (Am^ = m\ — m\ = (7.2 ■ ■ -8.9) x 10~ 5 eV 2 and |Am| 2 | = \m\ — 
m 2 \ = (1.7 ■■ -3.3) x 10~ 3 eV 2 ) and three neutrino mixing angles (30° < 9 12 < 38°, 
36° < 6 23 < 54° and 13 < 10°) at the 99% confidence level [13]. In addition, the latest 
cosmological constraint on the absolute values of m i is m 1 +m 2 + m 3 < 0.61 eV at the 
95% confidence level [14]. Since the unification of leptons and quarks is only possible 
at a superhigh energy scale (e.g., the GUT scale or the seesaw scale), it makes sense 
to calculate the running masses of three light neutrinos together with those of charged 
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fermions. We shall illustrate the running effects of m^fi) for // ^> M z by using the 
RGEs both in the SM and in the minimal supersymmetric standard model (MSSM). 

• The up-to-date values of quark masses [11] given by the Particle Data Group (PDG) 
at low energies are adopted, and a few new energy scales (such as /i ~ 1 TeV and 
10 12 GeV) are taken into account. Any models for new physics beyond the SM may 
introduce new energy scales above the Fermi scale, among which the fi ~ 1 TeV scale 
is most appealing because it will soon be probed by the LHG Possible new physics 
at this energy frontier is likely to be responsible for the origin of fermion masses and 
flavor mixing, or it can at least shed light on these fundamental problems. On the other 
hand, the intriguing seesaw [12] and leptogenesis [15] mechanisms have motivated a 
lot of neutrino physicists to study the properties of heavy Majorana neutrinos in order 
to simultaneously account for the lightness of three known neutrinos and the matter- 
antimatter asymmetry of our universe. The particularly interesting mass scales of such 
heavy Majorana particles are fi ~ 10 9 GeV to 10 12 GeV, where the flavor effects on 
leptogenesis [16] have to be treated carefully. Hence we shall make use of the RGEs to 
run the lepton and quark masses to these new energy scales, just for the convenience 
of model builders. 

• The treatment of quark masses crossing the flavor thresholds is improved. We shall 
use the matching conditions only at the flavor thresholds \i = mqimQ) (for Q = c,b, t) 
and calculate light quark masses at any other high energy scales with the help of the 
RGEs. This approach has been clearly described in Ref. [17]. In comparison, the 
running quark masses between two neighboring flavor thresholds were just computed 
with the matching conditions in Ref. [10]. The running- matching-running scheme 
proposed in Ref. [17] and used in our calculations is no doubt more reasonable. 

• The effective coupling constant a s , which governs the strength of strong interactions, 
is calculated by numerically solving its RGE. In contrast, the values of a s at different 
energy scales were evaluated via the analytical relation between a s and the asymptotic 
scale parameter A in Ref. [10]. The latter method may result in some unnecessary 
uncertainties due to the expansion of 1/ In (/i 2 /A 2 ), as shown in Refs. [17,18]. 

Our main numerical results will be tabulated, as done in Ref. [10], to serve for a useful 
reference in building specific models at various energy scales. 

The remaining parts of this paper are organized as follows. In Sec. II, we summarize the 
input data which include the current masses of three light quarks at fi = 2 GeV, the values 
of m c (m c ) and m b (m b ), the pole mass of the top quark, the strong gauge coupling a s and 
the fine structure constant a at M z . We also collect the relevant formulas for the RGEs and 
matching conditions of quark masses and a s . The strategy to deal with the flavor thresholds 
and evolve the charged fermion masses to the Fermi scale is outlined. Sec. Ill is devoted 
to the calculations of running fermion masses up to the GUT scale both in the SM and 
in the MSSM. Indeed, the true running parameters above the Fermi scale are the Yukawa 
couplings of leptons and quarks because these particles can only acquire their masses after 
the electroweak symmetry breaking (i.e., below the Fermi scale). The running masses of 
three light neutrinos are also illustrated for completeness. Finally, we make a brief summary 
of our main results in Sec. IV. 
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II. BELOW THE FERMI SCALE 



First of all, let us give some concise comments on the definition of fermion masses. It is 
important to specify the theoretical framework when discussing quark masses, since they are 
renormalization-scheme-dependent. There are two very common renormalization schemes. 
One of them is the on-shell scheme, in which the position of the pole is the definition of 
the physical mass M. The other is to define the running or renormalized mass m(n) in 
the dimensional regularization scheme with the modified minimal subtraction (MS), where 
H denotes the renormalization scale. Both definitions are suitable for charged leptons. But 
the situation is quite different in the quark sector due to the nonperturbative nature of 
the quantum chromodynamics (QCD) at low energies: light quarks are always confined in 
hadrons and can never be observed directly. Even for the heaviest quark, the top quark, it is 
impossible to completely eliminate the nonperturbative effect on the value of m t extracted 
from the direct measurements. Hence the pole mass of a quark is not well-defined. Note 
that the QCD Lagrangian has a chiral symmetry in the limit where all quark masses are 
vanishing. The scale of dynamical chiral symmetry breaking is A x ~ 1 GeV [19], which 
can be used to distinguish between light {m q < A x for q = u,d, s) and heavy (m q > A x for 
q — c,b, t) quarks. 

Another important point is the decoupling of heavy flavors. Because of the hierarchical 
mass spectrum of quarks, one should integrate out the heavy degrees of freedom when 
considering the properties of light flavors. Therefore, some consistent matching conditions 
should be established between the full and effective theories at the scales characterized by 
the masses of heavy flavors. We shall work in the framework elaborated on in Ref. [17]. 



A. Running Quark Masses 

With the help of the chiral symmetry, one may extract the values of m u /m d and m s /m d 
from the masses of pion and kaon in a way independent of the renormalization scale [20]. 
The mass of the strange quark can be determined from the spectral function sum rules or 
lattice QCD simulations [19,20]. The up-to-date values of m u , m d and m s given by the PDG 
are [11] 

m u (2 GeV) = 1.5 ~ 3.0 MeV , 

m d (2 GeV) = 3-7 MeV , (1) 
m s {2 GeV) = 95 ± 25 MeV . 

Note that these running masses are evaluated at /i — 2 GeV with three active quark flavors 
(u, d, s). The evolution of m q (/i) (for q — u,d, s, c, b, t) is governed by the following RGE: 

„.*^> = _ 7(a>>) = _g 7j ^y +1 m>) , (2) 

where a s (/x) = ^/(47r) is the effective coupling constant of strong interactions with g s being 
the strong gauge coupling. The values of r ) i (for % — 0, 1, 2, 3) are given by [21-23] 
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where n g is the number of active quark flavors with masses m q < /i] £(3) ~ 1.202057, 
((4) = vr 4 /90 « 1.082323 and C(5) « 1.036928 are the Riemann zeta functions. To find 
the solution to Eq. (2), we write out the detailed /z-dependence of a s (fj) in terms of the 
Callan-Symanzik beta function [11], 
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Note that we adopt the MS scheme for the RGEs throughout this paper. The solution to 
Eq. (2) can be expressed as [23,24] 



™ q {v) = K{a a (n))m q , 
where m q denotes the renormalization-invariant quark mass, and 



K(a s ) 
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In Eq. (7), the terms of O(aj) and smaller have been omitted and the coefficients C 123 are 
defined as 
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In a theory with multiple energy scales, such as the QCD with n l = Uq — 1 massless 
quarks and one heavy quark Q (for m Q ji), one should integrate out the heavy field and 
construct an effective theory by requiring its consistency with the full n^-flavor theory at 
the heavy quark threshold f/ n ^ = 0(m Q ). In this sense, Eq. (4) is valid between two quark 
thresholds, which are defined by = m q (m q ). Then we are in a position to address the 
decoupling of heavy quarks and consider the matching conditions at the flavor thresholds. 
In our calculations, we use the matching relation between the strong coupling constants of 
the neighboring flavors [17]: 



Oq- 1 )/ \ /-2 ( n q) I \ 
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where Q is already known at the three-loop level, i.e., 
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in the MS scheme. 

Now we consider the input values of three heavy quark masses. For charm and bottom 
quarks, their masses extracted from the heavy quark effective theory, lattice gauge theory 
and QCD sum rules are consistent with one another if they are all spelled out in the same 
scheme and at the same scale [11]: 



m c (m c ) = 1.25 ± 0.09 GeV , 
m b (m b ) = 4.20 ± 0.07 GeV . 



(11) 



The top quark can be directly measured because its lifetime is shorter than the typical time 
scale of nonperturbative strong interactions Aq^d [25], where the magnitude of Aq CD is 
about several hundred MeV. The pole mass of the top quark extracted from the average of 
several direct measurements is [11] 



M, = 172.5 ± 2.7 GeV 



(12) 
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Note that the nonperturbative contribution to the top quark mass may be proportional to 
Aq CD , which is much smaller than the present experimental error bar. The pole mass M q 
(for q — c,b, i) can be translated into the running mass at /i — M q through [24] 



M q = mJM. 
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.(3) 



(n ) 

of light quarks [24]. In addition to the matching condition of a s q in Eq. (9), we need to 
know the matching condition of the running quark masses at the flavor thresholds [17,26]: 
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where -B 4 ~ —1.762800. In practical calculations, we apply Eqs. (9) and (14) to the evolution 
of quark masses just at the thresholds fj,( n *\ Between two neighboring flavor thresholds, the 
RGE given in Eq. (2) or Eq. (6) will be used. We remark that this treatment is more 
reasonable than that adopted in Ref. [10]. 

It is worthwhile to mention that the running and decoupling of the strong coupling and 
quark masses have been built into the Mathematica package RunDec by Chetyrkin, Kiihn 
and Steinhauser [17]. We use the same formulas in our calculations and find that the results 
are in good agreement with those calculated by using RunDec, if the decoupling scale is 
chosen as /i (n <?) = m (m ). 



B. Running Charged-lepton Masses 

We proceed to discuss the running masses of charged leptons [27]. The /x-dependence of 
the fine structure constant a, including QCD corrections, is described by [28] 
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with Qj being the electric charge of a charged fermion (i.e., / = e,/i, r for charged leptons 
and / = u,d,s,c,b,t for quarks) and iV/ being the color factor (i.e., N l c — 1 for charged 
leptons and iV| = 3 for quarks), and 
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The pole masses of three charged leptons can be unambiguously measured in experiments 
and their values have been determined to an unprecedented degree of precision [11], 



M e = 0.510998918 ± 0.000000044 MeV , 
M„ = 105.6583692 ± 0.0000094 MeV , 

fj, 

M T = 1776.99l£;26 MeV . 



(18) 



One may convert the pole mass M t into the running mass m^/i) by using the following 
equation [29]: 



m^) = M { 1 - 
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(19) 



where the subscript / runs over e, ji and r, and the terms of 0(a 2 ) and smaller have been 
omitted. 



C. The Strategy 

The formulas given in the preceding subsections allow us to calculate the running masses 
of quarks and charged leptons up to the Fermi scale. The basic strategy of our numerical 
calculations is two-fold: 

1. We shall use the values of strong and electromagnetic coupling constants given by the 
PDG at M z [11], 

a s (M z ) = 0.1176 ±0.002 , 

a(Mz)- 1 = 127.918 ±0.018 , (20) 
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where M z = 91.1876 ± 0.0021 GeV is the mass of the Z boson and it lies in the range 
fj, > fi^ = m b {m b ). We shall only use the central value of M z in our calculations, 
because its error bar is negligibly small. From Eq. (4) with n q = 5 and Eq. (20), 
we obtain the corresponding strong coupling constant at the bottom quark threshold 

c4 5) (// 5) ) = 0.223to;oo7 witn A* (5) = m b( m b)- As for a ~ l > we can directly compute its 
values at some typical energy scales by using Eq. (15). Our numerical results, which 
will be used in the subsequent calculations, are listed in TABLE I. 

2. By using the RGEs with n q = 4, we are able to run three light quark masses and three 
charged- lepton masses to the first heavy flavor threshold / u ( - 5 * ) . Then we implement the 
matching conditions of a s and m q to cross this threshold. There is no flavor threshold 
on the way from fi^ to M z , and thus the relevant equations with n — 5 can be 
used for numerical calculations. We have also calculated the running masses mg(/i) of 
heavy quarks Q at the scale /i <C // n< ^ = rngfrng), which is below the flavor threshold 
^ n Q\ In our evaluation, the relevant equations with n Q active quark flavors have been 
used. For example, the running top quark mass at the Fermi scale m t (M z ) is actually 
computed by using the RGEs with n q = 6. 

Our numerical results for the running masses of quarks and charged leptons below the Fermi 
scale are given in TABLE II and TABLE III, respectively. We do not consider the running 
masses of three light neutrinos in this energy region, just because their interactions with 
other particles are too weak and their pole masses are too tiny. On the one hand, one 
expects that the changes of neutrino masses with respect to the energy scales are negligibly 
small. On the other hand, it is just the smallness of neutrino masses that hints at the possible 
existence of certain new physics scales (e.g., the seesaw scales /i ~ 10 9 GeV to 10 12 GeV). 
We shall illustrate the running effects of three neutrino masses from M z to fi ~ 10 12 GeV 
by using the RGEs of both the SM and the MSSM in the next section. 

III. ABOVE THE FERMI SCALE 

Above the Fermi scale, the unbroken electroweak gauge symmetry forbids leptons and 
quarks to acquire their masses. The actual meaning of a fermion mass in this energy 
region is a measure of the nontrivial Yukawa coupling eigenvalue y^. One commonly defines 
m f = yf v a b° ve the Fermi scale, just like the definition below the Fermi scale, where 
v pa 246 GeV is the vacuum expectation value of the neutral Higgs field in the SM. In the 
conventional seesaw models [12] 1 , the heavy Majorana neutrinos must be integrated out 
below the seesaw scale and the effective coupling matrix of three light Majorana neutrinos 
is given by the well-known seesaw relation k = —Y l/ M^ 1 Yj, where Y u denotes the neutrino 
Yukawa coupling matrix and M R is the right-handed Majorana neutrino mass matrix. The 



Here we follow the widely- accepted and well- motivated assumption that massive neutrinos are 
Majorana particles. The RGE running effects of Dirac neutrino masses and flavor mixing param- 
eters have been studied in Ref. [30]. 
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light neutrino masses are therefore given by m i = k^v 2 with n i (for % = 1,2,3) being the 
eigenvalues of k. Above the seesaw scale, the flavor threshold effects induced by the masses 
of heavy Majorana neutrinos have to be carefully treated and the details of Y v have to be 
model-dependently assumed [31]. Hence we shall only evaluate the running neutrino masses 
m i from the Fermi scale up to the seesaw scale, in order to avoid the complications and 
uncertainties associated with the seesaw thresholds. 

Then what we are concerned with is the evolution of four Yukawa coupling matrices Y u , 
Y d , Y[ and k, whose eigenvalues correspond to the masses of up-type quarks (m u , m c , m t ) , 
down-type quarks (m d ,m s ,m b ), charged leptons (m e ,m^,m T ) and neutrinos (m 1 ,m 2 ,m 3 ). 
Their one-loop RGEs in the SM and MSSM can be found in Refs. [30,31]. The two-loop 
beta functions of Y u , Y d and Y l have already been derived in Ref. [32]. Their lengthy 
expressions will not be quoted here, but they will be used in our numerical calculations. 
As for neutrinos, we only consider the one-loop RGEs for their running masses because 
our present knowledge on the absolute values of neutrino masses remain quite limited and 
uncertain [11]. The strategy of computing the running quark and lepton masses above 
the Fermi scale is outlined as follows. First, we use the quark masses and flavor mixing 
parameters obtained at the Fermi scale to reconstruct the Yukawa coupling matrices Y u and 
Y d . Second, the RGEs of Y u and Y d are solved and their eigenvalues are determined to give 
the running quark masses at every energy scale of our interest. We may follow a similar 
procedure to reconstruct Y l and k from current experimental data at M z , and then we solve 
their RGEs to find out the running lepton masses above M z . Because the RGEs of Y u , Y d , 
Y l and k are more or less entangled, our numerical calculations will be done simultaneously 
for quarks and leptons. 

Now that the running quark masses are to be evaluated in the supersymmetric framework, 
the matching procedure from the SM to the MSSM has to be taken into account. The RGEs 
in a supersymmetric theory are usually derived in the DR scheme based on the dimensional 
reduction and the minimal subtraction [33], while the experimental data on quark masses are 
extracted by using the MS scheme. Hence the transition between these two schemes needs 
to be treated in a consistent way [34,35]. To get around the occurrence of intermediate 
non-supersymmetric effective theories, here we follow Ref. [35] to adopt the common scale 
approach with all the supersymmetric particles being roughly at a common scale M. In our 
analysis, we set the decoupling scale to be \x = M = Mz- Furthermore, the input values 
of a s (M z ) and m JM Z ) in the MS scheme will be converted into those in the DR scheme 
at M z . First, since the two-loop RGEs will be used above the Fermi scale, we consider the 
transition of a s from the MS scheme to the DR scheme at the one- loop level [36], 



Given the input value of a s (M z ) in Eq. (20), a® R (M z ) ~ 0.1187 can be obtained and will 
be used as the input value in our numerical calculations. We take account of the matching 
between the SM and MSSM at M z ; i.e., a^ M = ( s a^ SSM . In the common scale approach, 
the decoupling constant approximates to £ s ~ 1 at the one-loop level, implying that the top 
quark and heavy sparticles are simultaneously decoupled. Second, the one-loop matching of 
quark masses between the MS and DR schemes is simply given by 




(21) 
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where a e is the evanescent coupling appearing in a non-supersymmetric theory renormalized 
in the DR scheme. We have a e ~ af R up to 0(a 2 s ) [35]. It is well known that the matching 
effects between the SM and MSSM are only significant for the masses of down-type quarks 
and charged leptons because of the large-tan (3 enhancement [37-40]. At the one-loop level, 
the threshold corrections read 

mssm _ m i (2 o\ 
m< "l + ^tan/T 

where % — d,s,b for down-type quarks or % — e, //, r for charged leptons, and the definitions of 
e i can be found in Refs. [39,40]. In our numerical analysis with the common scale approach, 
e i can be as large as one percent. Hence the supersymmetric threshold corrections are 
particularly relevant in the case of sizable tan/5 (e.g., tan/5 = 50). It is worth remarking 
that we have tried to avoid the details of any specific supersymmetric models within the 
scope of this work. If they are taken into account, however, a more careful treatment of the 
decoupling of supersymmetric partners will be unavoidable [37-39]. 

Before doing the numerical calculations, let us briefly summarize the relevant data to be 
input at the Fermi scale. 

1. In the basis where the Yukawa coupling matrix Y u is diagonal (i.e., Y u = Diag{y u , y c , y t } 
with y q = mjv for q — u,c,t), the mass eigenstates of down-type quarks (d, s, b) are 
related to their weak eigenstates (d', s', b') by the unitary Cabibbo-Kobayashi-Maskawa 
(CKM) matrix V [41]. Namely, V^Y d Y%V = Diag{y^,y^ y 2} with ^ = m j v for 
q = d,s and b. It is therefore possible to reconstruct Y d from its eigenvalues and V. 
Because quarks are Dirac particles, the CKM matrix V can be parametrized in terms 
of four independent quantities, such as the moduli of three independent elements of 
V and one CP- violating phase [42] . For example, one may choose the following set of 
parameters, which have been precisely measured, to parametrize V at M z [11]: 

\V cb \ = (41.6 ± 0.6) x 10" 3 , 
\V US \ = 0.2257 ±0.0021 , 
\VJ = (4.31 ± 0.30) xlO- 3 , 

sin 2j3 = 0.687 ±0.032 , (24) 



where (3 = arg [— {V^V* b ) / (V td V t * b )] is an inner angle of the CKM unitarity triangle. We 
shall only use the central values of the above parameters in our numerical calculations, 
because the running quark and lepton masses are actually insensitive to these inputs. 
On the other hand, the values of quark masses at the Fermi scale have been evaluated 
in Sec. II and listed in TABLE II. 



2. Without loss of generality, we choose the flavor basis where the Yukawa coupling 
matrix Y { is diagonal (i.e., Y l = Diag{y e , y , y T } with y l = mjv for I — e, /i, r). In this 
basis, the mass eigenstates of three light Majorana neutrinos (z/ 1? v 2 , u 3 ) are linked to 
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their weak eigenstates (u e , v , v T ) through the unitary Maki-Nakagawa-Sakata (MNS) 
matrix U [43] 2 . Namely, U^kU* = Diag{/t 1? k 2 , /%} with k { = m i /v 2 for i = 1,2 and 
3. It is then easy to reconstruct the symmetric matrix k from m i and U at the Fermi 
scale. A complete parametrization of the MNS matrix U needs three mixing angles 
and three CP- violating phases [45] , 

c 12 c 13 s 12 c 13 
U = I — S 12 C 2 3 — C 12 S23 S 13 e * 1 c 12 c 23 ~~ s 12 s 23 S 13 e * 1 S 23 C 13 \ \ ^ ® | (^5) 
S 12 S 23 _ C 12 C 23 S 13 e 1 ~ C 12 S 23 ~ S 12 C 23 S 13 e 1 

with Cj - = cos 6 i; j and s^- = sin^- (for ?j = 12, 13 and 23). A global analysis of current 
neutrino oscillation data yields the constraints [13] 

30° < 12 < 38° , 
36° < # 23 < 54° , 

13 < 10° (26) 




at the 99% confidence level. Three CP- violating phases ^, p and cr remain entirely 
unrestricted. It has been noticed that the running behaviors of neutrino masses are 
essentially insensitive to three neutrino mixing angles and three CP-violating phases 
in the SM or in the MSSM with small tan/3 [31]. Therefore, we shall simply assume 
5 l — p — a — and take the best-fit values 9 12 w 33.8°, 6 23 pa 45° and 9 13 pa 0° as the 
typical inputs at M z . Note that the absolute values of m i remain unknown, although 
their upper bound is expected to be 0{\) eV. For illustration, we only consider two 
possibilities in our numerical calculations: (A) the normal neutrino mass hierarchy 
with m 1 = 0.001 eV and m 1 < m 2 m 3 ; and (B) the nearly degenerate neutrino 
mass spectrum with m l = 0.2 eV and m 1 < m 2 < m 3 . The best-fit values Am^ = 
8.0 x 10" 5 eV 2 and Am 2 2 = 2.5 x 10~ 3 eV 2 are input at M z in both cases. This 
simplified treatment can easily be improved in the future, once more experimental 
data on m i and U are available. 

3. In the SM or MSSM with the gauge group S77(3) c x SU(2) L x U(1) Y , three gauge 
couplings g s , g and g' are given by 

gl = 47r «s , 

g 2 = 4vra/ sin 2 6 W , 

g' = gtiw9 w , (27) 

where 9 W is the weak mixing angle. In the SU(5) GUTs with or without supersym- 
metry, the gauge coupling constants g 1 , g 2 and g 3 are usually normalized as g 3 = g s , 



2 Because the seesaw scales are assumed to be extremely higher than the Fermi scale in this paper, 
U is expected to be unitary as an excellent approximation. See, e.g., Ref. [44] for a detailed analysis 
of the non- unitary neutrino mixing in the seesaw scenarios. 
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g 2 = g and g 1 = ^5/3 g'. The RGEs of g i (for i = 1,2,3) are given in Ref. [32]. The 
input parameters include [11] 

sin 2 e w (M z ) = 0.23122 ± 0.00015 , (28) 

as well as a s (M z ) and ct(M z ) given in Eq. (20). Besides the gauge couplings, the 
quartic Higgs coupling A = 2m 2 H / v 2 in the SM is also needed in our calculations. The 
RGE of A can be found in Ref. [32]. We shall typically take m H = 140 GeV for the 
Higgs mass, just for the sake of illustration. In the MSSM, the free parameter tan/3 
is defined as tan/3 = v u /v d with v d and v u being the vacuum expectation values of 
the Higgs fields which couple respectively to the down- and up-type quarks. Note 
that v 2 = v J + v 2 = AM^I g 2 ~ (246 GeV) 2 has been fixed by the measurements of 
the W-hoson mass and the electroweak gauge coupling. To be specific, we shall input 
tan/3 = 10 and tan/3 = 50 to illustrate the running fermion masses in the cases of 
small and large tan (3. 

With the help of the above inputs, we are then able to numerically solve the RGEs and 
obtain the running masses of quarks and leptons at various energies above the Fermi scale. 

Our numerical results at n = 1 TeV, 10 9 GeV, 10 12 GeV and A GUT ~2x 10 16 GeV 
are summarized in TABLE IV, TABLE V and TABLE VI, where the uncertainties of the 
output masses mainly come from those of the input masses at M z . Some brief comments 
on these results are in order. 

• In both the SM and the MSSM, a common feature of the outputs is that all the quark 
masses decrease with increasing energy scales. Nevertheless, their running behaviors 
are quite model-dependent. 

• TABLE IV and TABLE VI show that there exists a maximum for the running mass 
of each charged lepton in the SM. Our numerical analysis indicates that the maximal 
values of m e , and m T are all located around // ~ 10 6 GeV. In the MSSM, however, 
three charged-lepton masses smoothly decrease as the energy scale increases. 

• We have ignored large uncertainties of the input values in evaluating the running 
masses of three light neutrinos, just because the relevant experimental data are quite 
inaccurate and incomplete. Our numerical results, which rely on several assumptions 
made above, can only serve for illustration. We show the scale dependence of running 
neutrino masses in FIGs. 1 and 2 for two different mass spectra. Comparing between 
these two figures, we observe that the RGE effects on neutrino masses are more signifi- 
cant in the m 1 ^ m 2 ^ m 3 case than in the m 1 < m 2 <C m 3 case. Similar observations 
have been made by a number of authors before (see, e.g., Ref. [31]). 

Note that the running behaviors of other parameters, such as the CKM and MNS mixing 
angles, CP- violating phases and gauge coupling constants, can be obtained from our program 
in a straightforward way. But here we only concentrate on the outputs of running quark 
and lepton masses. 
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IV. SUMMARY 



In this paper, we have updated the running masses of quarks and leptons both in the SM 
and in the MSSM at various energy scales, including \x = M z , 1 TeV, 10 9 GeV, 10 12 GeV 
and A GUT . Our motivation is simple but meaningful: we want to provide a reliable and 
up-to-date table of the running quark and lepton masses for particle physicists. Such a table 
will be very useful for the analysis of hadronic physics at low energies and for the building 
of new physics models at superhigh energies. 

The differences between our work and the previous works (in particular, Ref. [10]) have 
been summarized in the introductory section of this paper. For instance, the central value 
of the strange quark mass m s at M z is about 55 MeV today, but it was about 93 MeV as 
given in Ref. [10] about a decade ago. Hence it makes sense to recalculate the values of m s 
and other fermion masses at different energy scales. In particular, we hope that our new 
results for the running masses of leptons and quarks may help the model builders to get new 
insight into the flavor dynamics at the energy frontier set by the LHC and in the exciting 
era of precision neutrino physics. 
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TABLES 



TABLE I. The gauge couplings a s (fi) and a(^) _1 at a few typical energy scales, where the first 
and second errors of a(/i) _1 come from the uncertainties associated with a{M z )~ x and m q (m q ), 
respectively. 





M K) = m q (m q ) (GeV) 


a s (fj,) 




4 


m c {m c ) = 1.25 


n QQ7+0.027 
u - 00 ' -0.024 


134.116 ± 0.018±g;ig? 


5 


m h (m b ) = 4.20 


n 99 o+0.008 


132.406 ± 0.0181^25 


6 


m t (m t ) = 163.6 


0.108 ±0.002 


127.073 ± 0.018 ± 0.023 



TABLE II. Running quark masses from = m c (m c ) to = m t (m t ), where we have taken 
the VL-boson mass to be M w = 80.403 GeV. The pole masses of three light quarks are not listed, 
simply because the perturbative QCD calculation is not reliable in that energy region. 





m u (fi) (MeV) 


m» (MeV) 


m» (MeV) 


mM (GeV) 


m b (fi) (GeV) 


m t (/x) (GeV) 


m c (m c ) 


o c; 7 +0.99 
z - ' -0.84 




1 1 1 + 31 


1.25 ±0.09 


x. qq+0.28 
°- yy -0.26 


384.8l^:| 


2 GeV 


9 9+O.8 
z - z -0.7 


5.0 ± 2.0 


95 ±25 


1-0718:11 


5 05+ - 16 


318.4l^:| 


m b (m b ) 


i -° b -0.60 


4.221};^ 


80 ±22 


o.9oi±g:ils 


4.20 ± 0.07 


259.81^ 


M w 


1 9Q+0-50 
i - zy -0.43 


2.93±J;lf 


56 ± 16 


626+ ' 084 

U.DZD_ 085 


2.92 ±0.09 


173.8 ±3.0 


M z 


i 9 7 +0.50 
l.Zl_ 0A2 


2.9011;?* 


55+ 16 
00 -15 


0.619 ±0.084 


2.89 ±0.09 


171.7 ±3.0 


m t (m t ) 


i 99+0.48 


2-7611:15 


52 ± 15 


0.590 ± 0.080 


2.75 ±0.09 


162.9 ±2.8 


M q 








1.77 ±0.14 




172.5 ±2.7 


m q (M q ) 








1 11+ - 11 


4.08 ± 0.08 


162.2 ±2.8 
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TABLE III. Running charged-lepton masses below m t (m t ), where the uncertainties of m z (//) 
are determined by those of M l . 





m e (n) (MeV) 


m M (/z) (MeV) 


m T (/x) (MeV) 


m c {m c ) 


0.495536319 ± +0.000000043 


104.47400561°;™ 


1774.90 ± 0.20 


m b (m b ) 


u.4yciuy4iyt>_ 000000043 


103.9951891 ± 0.0000093 


1767.08 ± 0.20 


M w 


0.486845675 ± 0.000000042 


102.7720886 ± 0.0000092 


1747.12±°; 2 ° 


M z 


0.486570161 ± 0.000000042 


102.7181359 ± 0.0000092 


1746.24 1°; 2 ° 


m t (m t ) 


0.485289396 ± 0.000000042 


102 4fi7si^^+ - 0000091 

1UZ.4D 1 olOO_ 000 Q 09 2 


1742.15 ±0.20 


M l 


0.510998918 ± 0.000000044 


105.6583692 ± 0.0000094 


1776.90 ± 0.20 
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TABLE IV. Running quark and lepton masses above M z in the SM with m H = 140 GeV, where 
the uncertainties of m^(fi) result from those of rrij(M z ). Here we have used A GUT = 2 x 10 16 GeV. 
Case A and case B represent two different neutrino mass patterns with m 1 (M z ) = 0.001 eV and 
m 1 (M z ) = 0.2 eV, respectively. 





H = M Z 


H = 1 TeV 


H = 10 9 GeV 


H = 10 12 GeV 


A* — A-GUT 


m u (ji) (MeV) 


-, 97 +0.50 
i - z ' -0.42 


i m +0 - 43 

i - iu -0.37 


fi7+ - 27 
u -°' -0.23 


58 +a24 
u.oo_ 20 


0.48+°$ 


m d (fi) (MeV) 


2.90^; 24 


O cn+1.08 
z - ou -1.03 


1-56183& 


1 Q4+0.60 




m» (MeV) 


55111 


4 7+l4 
4 ' -13 


30^ 


26±f 


22± 6 7 


m c (p) (GeV) 


0.619 ± 0.084 


n C09+0.074 
u.ooz_ 073 


n 097+0.048 
U - JZ ' -0.047 


o-28i±8:8g 


n 9 qc;+0.035 
U.ZOO_ 034 


m b (jj) (GeV) 


2.89 ± 0.09 


2.43 ± 0.08 


1.42 ± 0.06 


1.21 ±0.05 


1.00 ±0.04 


m t (fi) (GeV) 


171.7 ±3.0 


150.7 ±3.4 


99.llt° 


86.7±S:g 


74.0+ 4 ;? 


m e {n) (MeV) 


0.486570161 
±0.000000042 


0.495901601 
±0.000000043 


0.501014122 
±0.000000043 


0.490856087 

+0.000000042 
-0.000000043 


0.469652046 
±0.000000041 


m» (MeV) 


102.7181359 
±0.0000092 


104.6880645 

+0.0000094 
-0.0000093 


105.7673562 

+0.0000095 
-0.0000094 


103.6229311 

+0.0000092 
-0.0000093 


99.1466226 
±0.0000089 


m r (/x) (MeV) 


1746.24l°; 2 ° 


1779.74 ± 0.20 


1798.11+° JJ 


1761.67 ±0.20 


1685.58 ±0.19 


Case A 




( eV ) 


0.001 


0.001 


0.001 


0.001 




Am 2 ^) (eV 2 ) 


8.0 x 10" 5 


9.1 x 10~ 5 


1.3 x 10~ 4 


1.5 x 10" 4 




Am 2 2 (/x) (eV 2 ) 


2.5 x 10~ 3 


2.9 x 10~ 3 


4.2 x 10~ 3 


4.6 x 10" 3 


rs_< 


Case B 




miifj) ( eV ) 


0.20 


0.21 


0.26 


0.27 


r^j 


Am 2 ^) (eV 2 ) 


8.0 x 10" 5 


9.1 x 10~ 5 


1.3 x 10~ 4 


1.5 x 10~ 4 




Am^/i) (eV 2 ) 


2.5 x 10~ 3 


2.9 x 10 -3 


4.2 x 10" 3 


4.6 x 10~ 3 
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TABLE V. Running quark and lepton masses above M z in the MSSM with tan f3 = 10, where 
the matching effect between the MS and MSSM and the MS-to-DR transition effect on the input 
parameters at M z have been taken into account. 





H = M Z 


\i = 1 lev 


1 r»9 f ' \T 

/i = 10 GeV 


i nl2 /" 1 ir 

fi = 10 GeV 


A 4 ~~ GUT 


m u {fi) (MeV) 


1 ov+0.50 


r+0.45 
1 - io -0.38 


n ^+0.30 
°- 75 -0.25 


n co+0.26 

0.62I 021 


n /in+0.20 


( \ f~\/l \T\ 

m d {fi) (MeV J 


o nn+1.24 

2.90I 119 


o on+0.96 

2.20I 091 


1 O1+0.54 
L21 -0.51 


n nc+0.43 

0.96I 040 


n vn+0.31 
°- 70 -0.30 


m s(W (MeV) 


55_i5 


42 ± 12 


23 ± 7 


1 o+6 
18 -5 


13 ± 4 


m c (p) (GeV) 


0.619 ± 0.084 


„ ---+0.077 
u - 00 ' -0.076 


n oeq+0.053 
U.oDO_ 052 


n onQ+0- 046 

U.CiUCi_ 


A OQc+0.037 
U.Z^D_o.o36 


«w) ( GeV ) 


o o r\ i nn 

2.89 ± 0.09 


■ 1 ■ "i • > i n no 

2.23 ± 0.08 


1.30 ± 0.05 


1.05 ± 0.05 


r\ t-7f\ I A A ( 

0.79 ± 0.04 


m t {fi) (GeV) 


171.7 ±3.0 


161.0l| 6 7 


125.2^ 7 ; 5 


111.0±f;| 


92.2l?| 


m e (fi) (MeV) 


0.486570161 
±0.000000042 


0.418436115 
±0.000000036 


0.358332424 
±0.000000031 


0.327996884 

+0.000000028 
-0.000000029 


0.283755495 

+0.000000024 
-0.000000025 


m» (MeV) 


102.7181359 
±0.0000092 


88.3347018 
±0.0000079 


75.6468538 
±0.0000068 


69.2429377 
±0.0000062 


59.9033617 
±0.0000054 


m T (fi) (MeV) 


1746.24+°; 2 ° 


1502.25 ±0.17 


1288.68 ±0.15 


1180.38+g;^ 


1021.9510^2 


Case A 




m iW (ev) 


0.001 


0.001 


0.001 


0.001 




A 2 / \ / -\ 7-9 \ 

Am^(/x) (eV z ) 


8.0 x 10~ 5 


8.7 x 10~ 5 


1.0 x 10" 4 


1.0 x 10~ 4 




Ami(ii) (eV 2 ) 


2.5 x 10~ 3 


2.7 x 10~ 3 


3.3 x 10~ 3 


3.1 x 10~ 3 




Case B 




n»i(^) (eV) 


0.20 


0.21 


0.23 


0.22 




Am^/x) (eV 2 ) 


8.0 x 10~ 5 


9.1 x 10~ 5 


1.5 x 10~ 4 


1.6 x 10~ 4 




Ara^^) (eV 2 ) 


2.5 x 10~ 3 


2.7 x 10~ 3 


3.3 x 10~ 3 


3.2 x 10" 3 
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TABLE VI. Running quark and lepton masses above M z in the MSSM with tan (3 = 50, where 
the matching effect between the MS and MSSM and the MS-to-DR transition effect on the input 
parameters at M z have been taken into account. 





71 AT 

H = M Z 


H = 1 lev 


fi = 10 GeV 


1 nl2 f 1 \T 

H = 10 GeV 


A 

A* — GUT 


m u {fi) (MeV) 


1 0^+0.50 

1.27I 042 


i - io -0.38 


n 7C+0.31 
u - ' °-0.26 


n co+0.26 

0.62I 022 


n /iQ+0.21 


( \ f~\/l \T\ 

m d {fi) (MeV J 


2.90I 119 


i +0.67 
1 - Oi -0.63 


n qc+0.39 
°- 86 -0.36 


n cn+0.31 

0.69+ 029 


n c;i +0.23 
U - Oi -0.22 


m s (jjL) (MeV) 


55_i5 


on+9 

29-8 


16 ± 5 


13 ± 4 


10 ± 3 


m c (p) (GeV) 


A /"I A I /~V no A 

0.619 ± 0.084 


„ ^.p^+0.077 
u - 00 ' -0.076 


n oca +0.054 
°- 364 -0.053 


n on A +0.046 

0.304I 0045 


n 9Q7+0.037 
-0.036 


m b (ji) (GeV) 


2.89 ± 0.09 


1.54 ± 0.06 


0.96 ± 0.05 


n vn+0.05 

0.79I 004 


0.61 ± 0.04 


m t {fi) (GeV) 


171.7 ±3.0 


161.3 ±3.7 


127.0^; 4 


113.2.4+^ 


94.7™ 


m e (fi) (MeV) 


0.486570161 
±0.000000042 


0.286562894 
±0.000000025 


0.254747607 
±0.000000022 


0.235672689 

+0.000000021 
-0.000000020 


0.206036051 
±0.000000018 


m» (MeV) 


102.7181359 
±0.0000092 


60.4961413 
±0.0000054 


53.7834823 
±0.0000048 


49.7577528 

+0.0000045 
-0.0000044 


43.5020305 
±0.0000039 


m T (fi) (MeV) 


1746.24to: 2 g 


1032.61 ±0.12 


937.72 ±0.11 


875.31 ±0.11 


773.44 ±0.10 


Case A 




m iW (ev) 


0.001 


0.001 


0.001 


0.001 




A 2 / \ / -\ 7-9 \ 

Am^(/x) (eV z ) 


8.0 x 10~ 5 


8.7 x 10~ 5 


1.1 x 10" 4 


1.0 x 10~ 4 




Ami(ii) (eV 2 ) 


2.5 x 10~ 3 


2.7 x 10~ 3 


3.3 x 10~ 3 


3.2 x 10" 3 


r^j 


Case B 




"hM (eV) 


0.20 


0.21 


0.23 


0.23 


r^j 


Am 2 21 (p,) (eV 2 ) 


8.0 x 10~ 5 


1.7 x 10" 4 


7.1 x 10" 4 


8.3 x 10~ 4 




Ara^^) (eV 2 ) 


2.5 x 10" 3 


2.7 x 10" 3 


3.8 x 10~ 3 


4.1 x 10~ 3 
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FIGURES 



(a) MSSM, tanp=50 
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(b) MSSM, tanp=10 
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(c) SM, m=140GeV 
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FIG. 1. The scale dependence of three neutrino masses m 1 (solid line), m 2 (dashed line) and m 3 
(dotted line) in the case of a normal mass hierarchy with m 1 (M z ) = 0.001 eV and m 1 < m 2 <^ m 3 . 
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■ (a) MSSM, tanp=50 

0.28 - 



0.24 - 



0.20 




■ (b) MSSM, tan[3=10 



0.28 - 




(c) SM, m=140GeV 



0.28 - 
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FIG. 2. The scale dependence of three neutrino masses m l (solid line), m 2 (dashed line) and 
(dotted line) in the case of a nearly degenerate mass hierarchy with m 1 (M z ) = 0.2 eV and 
~ m 2 — m.3. 
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